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We have developed an accurate simulation model of the large 9 inch photomultiplier tubes (PMT) used in water- 
Cherenkov detectors of cosmic-ray induced extensive air-showers. This work was carried out as part of the devel- 
opment of the Offline simulation software for the Pierre Auger Observatory surface array but our findings may be 
relevant also for other astrophysics experiments that employ similar large PMTs. 

The implementation is realistic in terms of geometrical dimensions, optical processes at various surfaces, thin- 
film treatment of the photocathode, and photon reflections on the inner structure of the PMT. With the quantum 
efficiency obtained for this advanced model we have calibrated a much simpler and a more rudimentary model of 
the PMT which is more practical for massive simulation productions. We show that the quantum efficiency declared 
by manufactures of the PMTs is usually determined under conditions substantially different from those relevant for 
the particular experiment and thus requires careful (re)interpretation when applied to the experimental data or when 
used in simulations. In principle, the effective quantum efficiency could vary depending on the optical characteristics 
of individual events. 

Key words: photomultiplier tube, photocathode, Fresnel equations, thin-film, complex index of refraction, simulation 



1. Introduction 



The Pierre Auger Observatory is the largest detector 
^»built to detect cosmic rays with energies above 10 18 eV 



[1]. The surface detector part consists of over 1660 water- 



Cherenkov detectors distributed over an area of more 
than 3000 km 2 . While the longitudinal part of the shower 
development is measured with the fluorescence detec- 
tor, the muonic and electromagnetic lateral components 
are measured at ground level with the surface detector 
(SD). The individual SD stations consist of 12 1 of puri- 
fied water in a light-tight container with highly reflective 
and diffusive inner surface [1]. The Cherenkov light pro- 
duced by the penetrating charged particles is observed 
with PMTs. The signal levels from the PMTs are con- 
stantly calibrated against the average response to the at- 
J> mospheric background muons that trigger at the lowest 
threshold level [2]. Through dedicated experiments [2, 3] 
the response to atmospheric background muons arriving 
from all directions [4] has been related to the selection 
of atmospheric muons arriving predominantly vertically 
through the center of the SD station. All signals are thus 
measured in relative units of VEM (Vertical Equivalent 
Muon), i.e. relative to the average signal that would be 
produced by the vertical-centered muon. Although this 
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auto-calibration procedure successfully removes most of 
the systematics due to the detector changes, there are 
nevertheless some applications that require a certain de- 
gree of knowledge about absolute values. The absolute 
number of detected muons and the size of the electro- 
magnetic fraction in the signal are important quantities, 
e.g. in studies of hadronic interactions in air-shower cas- 
cades [5] or for the purposes of primary particle identi- 
fication [6, 7], just to name a few. We have developed 
a detailed simulation of the processes of light detection 
in the PMTs to reproduce known calibration data of the 
SD stations and allow comparison of real and simulated 
events. This study is part of the efforts to produce a 
complete SD simulation chain for the Pierre Auger Ob- 
servatory [8, 9, 10]. 

2. Photomultiplier tube 

The SD station has a cylindrical shape. Its height is 
1.2 m and its radius is 1.8 m. The Cherenkov radiation 
emitted in the water of the SD station is captured by 
three 9-inch PMTs floating on top of the water container, 
positioned 1.2 m away from the cylinder axis with 120° 
separation in azimuth. The PMTs have been produced 
by Photonis; the particular PMT model XP1805-PA1 used 
by the Pierre Auger Observatory differs from their stock 
model XP1805 [11] only in the additional output from the 
last dynode. Together with the usual anode output these 
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Figure 1: Schema of the advanced simulation model. The main com- 
ponents of the simulation set-up are the window glass, photocathode, 
mirror, effective absorption by the neck, and the multiplier structure. 

two are used for the monitoring of the dynode-anode 
ratio and the low and high gain signal acquisition [12]. 

The geometry of the PMT can be well approximated 
by an oblate spheroid with two equal semi-principal axes 
of 108.8 mm length and a third axis 79.4 mm long (outer 
dimensions). The thickness of the glass window is 2.5 mm. 
The glass is composed of borosilicate (80% SiC>2, 13% 
B2O2, and 7% Na 2 0) with an index of refraction of 1.47 
and a strong increase of absorption for wavelengths be- 
low ~ 300 nm. The photocathode consists of a thin ~ 
20 nm layer of bi-alkali metal (KCsSb) and has a wave- 
length-dependent complex index of refraction [13]. 

Based on this kind of geometry and material speci- 
fications of the PMT, we developed a simulation model 
that retains the basic geometry properties of the PMT 
while simulating in great detail the physical processes 
leading to the photoelectron emissions. 

3. Advanced simulation model 

The geometry set-up of the advanced simulation mod- 
el is shown in Fig. 1. The whole shape is approximated 
with a full oblate spheroid with a mirrored back wall 
reaching 30 mm below the center. The reflectivity of the 
mirror is set to 97% and only ideal specular reflection is 
implemented. As in the case of the real PMTs, the photo- 
cathode on the inner side of the glass is separated from 
the mirror coating by a transparent gap of 5 mm. 

3.1. Inner structure of the PMT 

The loss of photons in the extended neck leading to 
the PMT base is simulated by the 100% absorbing patch 
with a radius of 40 mm placed at the top of the PMT. 
The multiplier structure (dynode stack) reaches well into 
the center of the PMT and is in this particular PMT en- 
closed in a metallic shield case of cylindrical shape. In 
the photon-tracking experiments we have found that the 
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Figure 2: Real (full line) and negative imaginary (dashed line) parts of 
the complex refraction index of the photocathode as a function of the 
wavelength. The triangles are experimental data points extracted from 
[13]. The crosses are extrapolated values. The dotted line represents 
the Fermi function used to model the transmission factor of the PMT 
glass. 

proper inclusion of this multiplier volume can greatly in- 
fluence the number of reflected photons (see Fig. 4 for a 
clear demonstration on how the parallel beam of light 
gets focused into the multiplier volume). The multi- 
plier structure is thus simulated with a cylinder of height 
60 mm and radius 16 mm. It is made of a copper-like 
metal and its effective reflectivity is set to 30% (for sim- 
plicity specular reflection only). 

3.2. Optical properties of the PMT glass 

The material of the glass window limits the spectral 
sensitivity in the short wavelength region. The borosili- 
cate glass has a cut-off wavelength of 270 nm (the point 
where it decreases below 10%). We have modeled this 
property with a simple transmission factor Tgi ass that is 
imposed upon the entering photons, 

T g ia S s = exp[-d/L(A)] = F(A; A m , AA), (1) 

where d is the thickness and L(A) is the absorption length 
of the PMT glass. The final transmission factor is mod- 
eled with a Fermi-function dependence on the wavelength, 
F(A; A t , AA) = 1/ (exp[(A t — A)/AA] +1) with A t = 285 nm 
for the transition wavelength and AA = 8 nm for the tran- 
sition width [14] (see dotted line in Fig. 2). 

3.3. Photocathode as a thin film 

The sensitivity of the different types of photocathodes 
is restricted by the photo-emission threshold for the long 
wavelengths and can vary with thickness at the short 
wavelengths. The bi-alkali component KCsSb is practi- 
cally the universal choice of photocathode material for 
Cherenkov light applications, and although it has been 
known for almost five decades, the availability of reli- 
able measurements of index of refraction has been rather 



scarce. In our simulation of the angular and wavelength 
dependencies of the quantum efficiency we have used 
the data points from the experimental compilation of the 
complex index of refraction from [13], where the optical 
properties of this bi-alkali material of similar thickness 
have been reviewed for the purposes of the solar neu- 
trino detector SNO. Due to the limited number of dif- 
ferent wavelength measurements we have used a simple 
linear interpolation of the complex index of refraction 
with wavelength extension along the lines of newer ex- 
perimental results from [15] (see Fig. 2). 

3.4. Photoelectron production 

When the photon reaches the thin layer of the photo- 
cathode it can either get reflected with probability R, it 
can get transmitted across the layer with effective proba- 
bility T, or it can get absorbed by the photocathode mate- 
rial with probability A = 1 — R — T. If the photon energy 
in the latter case is larger than the needed exit work, the 
electron can leave the photocathode. The corresponding 
quantum efficiency q pc of the photocathode surface can 
thus be written as 



<7pc(A,0) = A(A,0); 



(2) 



where A (A, 6), the absorption coefficient, depends on in- 
cidence angle 6 and wavelength A, and p C onv is the con- 
version factor. The conversion factor p conv does not de- 
pend on the angle of incidence [16, 17] but in principle it 
can still be wavelength dependent [18]. Due to the rela- 
tively narrow window of relevant wavelengths imposed 
by the water absorption (see Fig. 3) we approximate the 
conversion factor p con v for the purposes of our simula- 
tions with a constant. In section 3.7, we determine its ef- 
fective value from the separate simulation of the quantum 
efficiency experiment, along the lines of the experiments 
usually performed by the PMT manufacturers. In the 
next section we will derive the second missing parame- 
ter A(A,0). 

At this point it is worth to mention that the quantum 
efficiency q pc from Eq. (2) is the quantum efficiency of the 
photocathode surface and is only indirectly related to the 
overall quantum efficiency of the PMT, as found in the 
specifications of the manufacturers. The main difference 
comes from the fact that PMTs are designed to increase 
light collection through reflections on the mirror back 
wall. This increase is to some degree contained in the 
PMT specification. Nevertheless, its magnitude depends 
substantially on the particular way the quoted efficiency 
is obtained experimentally. Therefore, the quoted effi- 
ciency specifications of the PMT always require proper 
interpretation when used for a photocathode simulation. 
This fact is in general relevant also for other physics ex- 
periments involving PMTs. 

3.5. Thin-film Fresnel equations 

Studies on PMTs usually follow [19] to derive the an- 
gular and wavelength dependence of the optical proper- 



ties. Here, we follow the more concise derivation from 
[20] that is specifically dedicated to thin film treatment. 
Furthermore, we adopt from [20] the sign convention for 
the imaginary part of the index of refraction (see Fig. 2). 

Although our PMT has only one layer of thin-film 
photocathode, we will for the sake of clarity derive Fres- 
nel equations for a stack of m thin-film surfaces with cor- 
responding thicknesses di and indexes of reflection ny. 
The indexes n; vary with different wavelengths but in 
the expressions below we do not explicitly notate the de- 
pendence on A for the sake of clarity. This stack of thin 
films is on the in-going and out-going side surrounded 
by two semi-infinite substrates with indexes Mi n and n ut/ 
respectively. 

The angles of refraction in the consequent layers can 
all be obtained from the incidence angle 8[ n from Snell's 
law 

(3) 



n j cos Uj = 



n\ n sin 2 e^. 



This expression should be used for all m layers j — 1 . . . m 
even when it produces a complex cosine both in case 
of a complex index ft; or a total internal reflection n; < 



, sin #i n . Phase change in the layer j is denoted by 

d: 

Si = 2n-^ni cos 9;. 
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In the case of thin layers, the transmitted and reflected 
light combine coherently and the vector of normalized 
emergent fields f ou t can be expressed as in [20] with char- 
acteristic matrices of m stacked layers as 
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where /g and fu are normalized electrical and magnetic 
fields, respectively, and m, t] out are the tilted optical ad- 
mittances. Using short-hand notation for scalar products 
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the reflectance is expressed as 
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and the transmittance as 
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Absorptance follows from 
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A = An h 



(8) 



(9) 



The three quantities obey the relation R + T + A = 1. 
Note that while rj m is always real, rjj can be complex due 
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Figure 3: Wavelength dependency of the water absorption length L and 
the reflectivity R of the container walls. Data taken from [21] and [22], 
respectively. 

to absorption, and rji, t] out can become complex when 
the incidence angle 6- m increases above the angle of total 
internal reflection 0' ot defined by sinf^ ot = rij/n m . 

The equations above have to be considered for two 
polarization cases. Using the convention where p-polariz- 
ation implies a magnetic field component parallel to the 
interface boundary (TM wave) and s-polarization implies 
a parallel electric component (TE wave), two variants of 
the admittance emerge, 

rjj = tij/ cos dj and rf- = rij cos 9j, (10) 

where cosfy is again obtained from Eq. (3). The same 
polarization cases should also be used on rj^ and n oui . 
The upper expressions in (7), (8), and (9) thus have to be 
separately evaluated for the p- and s-polarization; hence 
the obtained reflectance, transmittance, and absorptance 
change accordingly into Rp /S , Tp /S , and A ps . For un- 
polarized light we can define average quantities R = 
l(R p +R s ),T= i(T p + T s ),andA = \{A V + A s ) which, 
as before, satisfy R + T + A = 1 . 

Cherenkov light from the injected muon is polarized 
perpendicularly to the emission cone. The rotational sym- 
metry around the muon track and the fact that after sev- 
eral reflections on the inner walls the photon polariza- 
tion is after ~ 30 ns completely randomized, enable us to 
use the above polarization-averaged expressions for re- 
flectance, transmittance, and absorptance. 

3.6. Simulation set-up 

The simulation of the PMT has been implemented 
within the Auger Off line software framework [23] where 
the water-Cherenkov detector simulation is implemented 
with the Geant4 toolkit [24]. The production of the 
Cherenkov light by the injected vertical-centered muons, 
the consequent absorption of the photons in the water, 
and the tracking of the photons and their reflections on 




Figure 4: Simulation of the Photonis quantum efficiency experiment for 
estimation of the conversion factor. Absolutely calibrated parallel beam 
of light is sent to the front of the PMT (photocathode is at bottom of the 
oval). Due to the reflection on the ellipsoidal mirror at the back of the 
PMT, the light gets approximately focused at the position of the mul- 
tiplier (caustics at the sides of the multiplier). The inclusion /exclusion 
of the multiplier volume can have great effect on the amount of the 
light returned to the photocathode surface and thus has to be properly 
simulated. In this figure the multiplier (c.f. white patch at top cen- 
ter) is obscuring most of the returning light in this configuration of the 
incoming beam. 

the container walls have been assigned to the Geant4 
part of the simulation with parametrizations shown in 
Fig. 3. The entry of the photons into the glass of the PMT 
window was also handled by Geant4. From this point 
on, the photons were tracked by our simulation model, 
based on the considerations made in the previous sec- 
tions. Using the index of refraction parametrization from 
Fig. 2 for n\, the individual photons were reflected from 
the photocathode with the probability following from the 
(polarization-averaged) reflectance in Eq. (7), or transmit- 
ted according to Eq. (8) into the vacuum (n ut = 1) of the 
PMT. The simulation proceeded with the photon track- 
ing inside the PMT which can, upon successful reflec- 
tions on the inner structure, produce another crossing of 
the photocathode (this time with the exchanged values 
for the "in" and "out" variables in the upper equations) 
and reentry into the water. In case of an absorption of the 
photon upon the crossing of the photocathode, the asso- 
ciated photoelectron is produced according to the quan- 
tum efficiency in Eq. (2). 

3.7. Calibration of the conversion factor 

For the wavelength of 375 nm, Photonis [11] quotes 
a value of quantum efficiency of the PMT as ^p m t 
0.26. They perform the measurement with an absolutely- 
calibrated parallel beam of light with a diameter of 170 mm 
illuminating the front of the PMT along the PMT axis. 
Correcting for the 70% collection efficiency (obtained from 
their own electric field simulations [14]), they consequently 
convert the measured current into a released photoelec- 
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Figure 5: Normalized distribution of photons arriving to the photo- 
cathode as a function of cos 28 where 8 is the incidence angle, for two 
populations, incoming (i.e. entering the PMT, full line) and outgoing 
(i.e. exiting the PMT, dashed line) photons. Due to the increased re- 
flectivity of the water-glass interface a depletion from the isotropic dis- 
tribution is observed for large incidence angles (cos(2 x 90°) = —1) 
of incoming photons. For the outgoing photons the small incidence 
angles (cos(2 x 0°) = 1) are suppressed due to the reflection and ab- 
sorption on the inner structure of the PMT. Note that the distribution 
for isotropic arrivals as a function of cos 28 is constant and equals to 
1/2. 



tron flux O pe . Their definition of PMT quantum effi- 
ciency is then ^p m t = Ope / <3? where <J> is the incident 
photon flux from the absolutely calibrated beam. 

We have reproduced this exact experimental setup 
and the beam geometry in our PMT simulation (see Fig. 4) 
in order to "reverse engineer" the PMT quantum effi- 
ciency <7pmt into the photocathode quantum efficiency 
<7p C . While using the surface quantum efficiency expres- 
sion from Eq. (2) we simulated a large number N of the 
375 nm photons in the incident beam and recorded the 
number of released photoelectrons N pe . To reproduce 
the quoted PMT quantum efficiency ^p m t = 0.26 with 
our simulated fraction Np e /N, the conversion factor in 
Eq. (2) had to be set to p C onv = 0.405. Thus, in effect, at 
A = 375 nm slightly more than 40% of the absorbed pho- 
tons in this geometry get converted to photoelectrons. 

3.8. Quantum efficiency 

Using the obtained conversion factor p CO nv = 0.405, 
the PMT simulation is run for a large number of in- 
jected vertical-centered muons. Fig. 5 shows the distri- 
bution of the incidence angle on the photocathode for 
photons entering (incoming) and exiting (outgoing) the 
PMT. On each crossing of the photocathode, absorbtance 
A (A, 9) is obtained form the Eq. (9) and the photoelec- 
tron is released with probability A(A,6) p C onv, as given 
by Eq. (2). We can derive the overall, incidence-angle av- 
eraged quantum efficiency of our PMT simulation model 
by dividing the number of released photoelectrons N pe 
by the number of photons N reaching the photocathode 



Figure 6: Quantum efficiency q (full line) defined as the probability 
for a photon to free an electron at the photocathode crossing. The 
dashed line corresponds to the effective quantum efficiency q e ff, i.e. 
the probability for an incoming photon to eventually free an electron. 
The effective quantum efficiency is used in the implementation of the 
simplified PMT model. 



(from any side) for the different wavelengths, 



Npe(A) 
N(A) • 



(11) 



The resulting quantum efficiency is shown in Fig. 6 (full 
line). This quantum efficiency, obtained indirectly from 
the refractive index parametrization in Fig. 2, reproduces 
well the known features of other experimental data [15] 
and the specification of the PMT manufacturer [11]. Fur- 
thermore, it gives the correct estimates of the average 
number of photoelectrons released by the traversing muon 
as measured in scintillator-triggered experiments with 
Auger SD stations [2, 3]. 

4. Relevance to other astrophysics experiments 

There are several important points that can be ex- 
tracted from the previous sections. For optical inter- 
faces, the Geant4 framework implements only the fa- 
miliar Fresnel equations of geometrical optics and does 
not include the physical optics relevant for a thin-film 
photocathode, as given by Eqs. (5-9). The dependence 
of the photocathode absorption (and consequently of the 
quantum efficiency) on the incidence angle is correctly 
described only by the thin-film equations which have 
to be implemented independently of Geant4. The ab- 
sorption increases up to 50% for incidence angles larger 
than the angle of total internal reflection for the glass- 
vacuum interface. No such increase is observed for the 
reverse, vacuum-glass transition of the photons. Further- 
more, the quantum efficiency obtained from the manu- 
facturer's calibration setup already contains a nontrivial 
fraction of the efficiency increase due to the reflectivity 
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Figure 7: Photon spectrum reaching the photocathode for 10000 in- 
jected vertical-centered muons. The main features of the curve are 
due to the shape of the Cherenkov emission spectrum, the wavelength 
dependency of the water absorption length (hump at ~ 420 nm), the 
wavelength dependence of the reflectivity of the inner walls, and the 
PMT glass transparency (cut-off at wavelengths < 280 nm). 



of the inner structure of the PMT. The "intrinsic" quan- 
tum efficiency of the photocathode layer thus has to be 
reverse-engineered by a procedure similar to the one de- 
scribed in Section 3.7. Clearly both effects are most rel- 
evant in cases of experiments with a distribution of in- 
coming photons substantially different from those char- 
acteristic of manufacturer's experimental setup. 

This may be of relevance to experiments employing 
similarly large PMTs in various Cherenkov detectors, like 
Amanda, IceCube, and IceTop [25, 26, 27], Kamiokande 
[28], NESTOR, NEMO, and Antares [29, 30, 31], Mini- 
BooNE [32], Xenon [33], Baikal and Tunka [34], and north- 
ern Auger Observatory [35], just to name a few. 

5. Simplified simulation model 

In a simulation, a vertical-centered muon injected into 
the water of the SD station has a 1.2 m tracklength and 
typically releases ~ 67500 Cherenkov photons in the en- 
ergy range between 1.5 and 5eV (250 to 830 nm in wave- 
length; see Fig. 7 for spectrum). Additional ~ 10% of 
photons get produced by the secondary delta rays, in to- 
tal giving <~ 73 100 photons per muon. All these photons 
require a complete tracking inside the water container by 
Geant4. The photons on their way undergo the absorp- 
tion in the water medium, get reflected and absorbed by 
the container walls. On average, only ~ 1340 (1.8%) of 
all photons eventually reach the PMT photocathode. As 
described in the previous sections, the photon tracks are 
then simulated inside the PMT and the probability of a 
photoelectron release is evaluated. A large fraction of 
produced photons will never reach the PMT and addi- 
tionally a fraction of the incoming photons will be ab- 
sorbed without producing any photoelectrons. To avoid 
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Figure 8: Ratio of outgoing vs. incoming photons as counted at the 
PMT photocathode. 



the time-consuming simulation of these photons that do 
not contribute to the PMT signal, a simplified simulation 
has been devised. The goal is to estimate the effective 
probability for an emerging photon to contribute to the 
signal and eliminate idle photons from the simulation. 
The corresponding statistical thinning can thus be applied 
already at the time of their production. 

5.1. Effective quantum efficiency for incoming photons 

After several reflections on the walls, the swarm of 
photons behaves similarly to a photon gas. As can be al- 
ready seen from the incidence-angle distribution of pho- 
tons on the photocathode (see Fig. 5), the distribution of 
arrival directions is close to isotropic. Instead of the de- 
tailed simulation of photon tracks inside the advanced 
PMT model we can reduce the complexity of the pro- 
cedure by introducing the simplified model of the PMT. 
The quantum efficiency in Eq. (11) stemming from Fres- 
nel equations of the advanced PMT model is replaced 
with an effective quantum efficiency q e g that describes the 
probability for an incoming photon to produce a photo- 
electron, 



</eff( A ) 



N pe (A) 



(12) 



Win (A) ' 

The ^ e ff is obtained from the advanced model by nor- 
malizing the statistics of the released photoelectrons by 
the number of incoming photons. The actual wavelength 
dependence of q e ft (A) is shown in Fig. 6 (dashed line). 

In the advanced model only 4.5% of the incoming 
photons return to the water medium, i.e. only 0.08% of 
all produced photons. In Fig. 8, a wavelength-resolved 
fraction of outgoing vs. incoming photons is shown. The 
fraction is as low as 2% for short wavelengths and does 
not exceed 22% for long wavelengths. 

Based on these facts we have developed an imple- 
mentation of the simplified PMT model that uses only 
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Figure 9: Representation of the simplified model of PMT. 
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Figure 10: Correction factor q e ff (A)/g(A) as a function of the wave- 
length. For short wavelengths the correction is of the order of 13% and 
goes up to 35% for long wavelengths. 



a rudimentary geometry as shown in Fig. 9. In this sim- 
plified model the photons are detected on contact with 
the outer surface and are consequentially removed from 
the simulation. No tracking of the photons on the in- 
ner structure of the PMT is performed. The chance of 
producing a photoelectron on the second crossing of the 
photocathode in the advanced model is accounted for in 
the simplified model by the usage of the effective quan- 
tum efficiency from Eq. (12). The relative ratio of the 
quantum efficiencies from Eqs. (12) and (11) is shown in 
Fig. 10. The ratio is mostly close to ~ 10% increase but 
then gradually climbs to ~ 35% for long wavelengths. 
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Figure 11: Comparison of the time dependence of the muon signal in 
an SD station for the advanced PMT model (full line) and the simplified 
version (dashed line). The main feature is an approximate exponential 
time dependence with a decay time T ~ 65 ns. The two models agree 
well with respect to the small details of the early signal and the latter 
decay part. A small discrepancy, limited to less than 5%, is observed 
only in the late tail. Note that the oscillating structure for time < 30 ns 
is due to the nonuniform arrival of photons after the first few reflections 
when photons can not be treated yet as homogeneous and isotropic 
photon gas. 



account also the losses due to the PMT collection effi- 
ciency this number increases to ~ 6.5. 

In Fig. 11 the results for the muon response in SD sta- 
tion is compared for both models. The main feature is an 
exponential decay but with a slightly changing exponent 
due to the wavelength dependence of the water absorp- 
tion and wall reflectivity. The characteristic decay time 
of the muon signal is around t w 65 ns. Except in the tail 
of relatively large times, t » t, where the discrepancy 
is limited to within 5%, the simplified model reproduces 
well the details of the muon signal from the advanced 
PMT model, especially the overall exponential decay and 
an oscillatory behavior for t < 20 ns. It is worth mention- 
ing here that the signal from the PMT in an SD station is 
sampled with a 40 MHz (25 ns) FADC, i.e. all these details 
will lie within one sampling bin. 



5.2. Statistical thinning 

Similar to the well established method in extensive 
air-shower simulations [36, 37], the concept of effective 
quantum efficiency enables us to implement efficient sta- 
tistical thinning of the simulated photons. At the mo- 
ment of production, out of N pr (A) photons with wave- 
length A only ^ e ff (A)Np r (A) photons with weight l/q e {{(\) 
undergo further simulation. As a result, such thinned 
photons upon reaching the photocathode have to pro- 
duce a photoelectron with probability 1. The average 
quantum efficiency (in photon energy scale where Che- 
renkov spectrum is flat) is q e ff — 0.22, resulting in a sim- 
ulation speed-up for almost a factor ~ 5. If we take into 



6. Summary 

We have implemented an advanced model for the 
simulation of the PMT response in a water-Cherenkov 
station of the surface detector. The model is based on the 
correct thin-film treatment of the photocathode optical 
processes with complex index of refraction. The model 
also includes tracking of the photons on the inner mirror 
surfaces and multiplier structure. 

To relate the obtained quantum efficiency of the pho- 
tocathode to quoted efficiency from the PMT specifica- 
tion we have in simulation reproduced the experimental 
set-up of the measurement done by the manufacturer. 
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This properly accounts for the increase in photon collect- 
ing due to the reflections on mirrored surfaces and inner 
structure of the PMT, and gives the average value of the 
conversion factor. 

Since a relatively large number of photons is involved 
in the creation of the muon signal, certain details of the 
simulation, like the dependence of the absorption on the 
incidence angles and the exact tracking of the photons on 
the inner structure, are washed out and can be treated in 
a phenomenological way. 

To reduce the simulation time as well as the complex- 
ity of the problem we have implemented another, sim- 
plified model of the PMT. Results for the average pho- 
toelectron release probabilities for different wavelengths 
in the advanced PMT model were in turn used to cali- 
brate the effective quantum efficiency of the simplified 
model. Results from both models agree well with the 
experimentally obtained number of produced photoelec- 
trons per muon injection. Finally, we have shown that 
the time dependence of the muon signal is almost not 
affected by the simplification in the simulation strategy. 
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